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Abstract 



We have investigated the magnetoplasmon excitations in arrays of cir- 
cular and noncircular quantum dots within the Thomas-Fermi-Dirac-von 
Weizsacker approximation. Deviations from the ideal collective excitations 
of isolated parabolically confined electrons arise from local perturbations of 
the confining potential as well as interdot Coulomb interactions. The latter 
are unimportant unless the interdot separations are of the order of the size of 
the dots. Local perturbations such as radial anharmonicity and noncircular 
symmetry lead to clear signatures of the violation of the generalized Kohn 
theorem. In particular, the reduction of the local symmetry from 50(2) to 
C4 results in a resonant coupling of different modes and an observable anti- 
crossing behaviour in the power absorption spectrum. Our results are in good 
agreement with recent far- infrared (FIR) transmission experiments. 
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I. INTRODUCTION 



Advances in nanofabrication technology have opened up the possibihty of studying elec- 
tronic systems under a wide range of confinement conditions. The systems of interest are 
usually fabricated by starting with a two-dimensional electron gas and applying various lat- 
eral modulation techniques.uH In this way, systems such as wires, dots, rings and antidot 
arrays have been realized. A powerful method for studying these systems is FIR absorption 
spectroscopy. The observation of collective plasmon excitations and their dispersion in an 
applied magnetic field can yield a great deal of information about the effects of electron- 
electron interactions and the influence of different forms of geometrical confinement. 

In this paper we are primarily interested in the magnetoplasmon excitations in quantum 
dot arrays. In particular, we focus on the experimental results of Demel et a/.@'i which reveal 
interesting features in the absorption spectrum as a function of magnetic field. It is by now 
well understood that the observed deviations from the ideal spectrum expected for parabolic 
confinementihl are mainly due to anharmonic perturbations of the confining potential. 
For parabolic confinement, the only allowed dipole excitations correspond to center of mass 



motion and have the frequencies Uj- = + {ujc/2y ± uJc/2, where Uq is the harmonic 
confining frequency and Uc is the cyclotron frequency. In a previous papeJll it was shown 
that an axially symmetric perturbation leads to a coupling of the center-of-mass mode to 
other dipole excitations of the harmonic potential, and gives rise to a satellite peak which 
tracks along the u!+ excitation. The calculated absorption spectrum was found to be in good 
agreement with experiment with regard to both the position and strength of the satellite 
peak. 

One feature which this earlier calculationlli did not account for was an anticrossing be- 
haviour observed at lower magnetic fields. It was originally conjecturedil that this feature 
was due to the nonlocal dynamic response associated with the finite compressibility of an 
electron gas and responsible for the wavevector dispersion of plasmons in a uniform sys- 
tem. An explanation on the basis of a noncircular confining potential was dismissed on the 
grounds that the observed anticrossing did not conform with its expected dependence on 
A^, the number of electrons in the dot. Subsequent work, however, shifted towards the non- 
circular confining potential as being the most likely explanation of the observed behaviour. 
Gudmundsson and Gerhardtsi performed random-phase approximation (RPA) calculations 
for a parabolic dot with 10 < A^ < 30 and argued, on the basis of the coupling between 
the center of mass and relative degrees of freedom induced by a x^y^ perturbation, that 
an anticrossing of the dipole mode with modes of higher multipolarity is to be expected. 
Pfannkuche and Gerhardts0 arrived at similar conclusions based on an exact diagonaliza- 
tion of the dot Hamiltonian with quartic perturbations, but only for A^ = 2. A classical 
electrostatic model was used by Nazin et aZ-IHI to study the magnetoplasmon modes in a 
parabolic dot with a cos(46') perturbation. This approach, which treats the electrons as 
a charged classical fluid, can be formulated in terms of an approximate energy functional 
which includes only the effects of the confining potential and the electrostatic repulsion 
between the electrons. The model is useful for large A^ and was shown to provide a good de- 
scription of the dynamics in anharmonic dots.0 However, the theory proved to be ill-defined 
for noncircular dots as a result of the singular behaviour of the equilibrium density at the 
edge of the dot,0 and a rigorous analysis of the dynamics could not be provided. 
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Our purpose here is to re-examine circular and noncircular quantum dots within the 
Thomas-Fermi-Dirac-von Weizsacker (TFDW)hydrodynamic theory developed in the con- 
text of parabolically confined electron slabs.EJ The theory was later applied to electron 
rme and was shown to give a good account of the edge magnetoplasmon modes observed 
in this system. The theory improves on the classical electrostatic model in that the quantum 
kinetic energy is included at the level of the TFDW approximation. In addition, exchange- 
correlation effects can also be included. Unlike other hydrodynamic-based models,EHli the 
TFDW hydrodynamic theory incorporates the calculated equilibrium properties of the elec- 
tronic density and thereby provides a more physical description of the collective modes. 

The approach that we use is best suited to treating the FIR response of dots with a 
large number of electrons where the response is expected to be dominated by collective-like 
excitations. For few electrons, the underlying electronic structure can play an important role. 
In this situation, methods based on the random phase approximation are more appropriate, 
and calculations for noncircular dots in this regime have recently appeared.0@ However, 
these calculations cannot be extended readily to dots with a large number of electrons which 
is the regime of interest in this paper. 

Our paper is organized as follows. In Sec. |I| we present an overview of the TFDW 
formalism and develop the constitutive equations for the theory in the context of a bi- 
directionally modulated electronic system. In Sec. we apply our theory to dot arrays and 
examine the effects of confinement symmetry and interdot coupling on the magnetoplasmon 
modes of the system. Finally, in Sec. we present our concluding remarks. A brief 
account of this work will appear as a conference proceeding.!! 



II. THE TFDW FORMALISM 

The equilibrium and hydrodynamic properties of a periodically modulated 2DEG within 
the TFDW approximation have already been discussed at some length in our previous work 
on uniaxially modulated electronic systems.c2l In this section, we provide an overview of 
the TFDW formalism and extend our earlier work to the more general case of bi-directional 
modulation. Since these calculations parallel those in Ref. [ ^], we only present the essential 
ingredients of the model here and refer to our earlier papers for details. 

A. Equilibrium Properties 

As in all density functional theory (DFT) schemes,!! the equilibrium properties of the 
system are obtained by finding the variational minimum of an energy functional. In atomic 
units (e^/e = m* = h = 1), the TFDW functional is given by 



E\n]= I dh 



n 2 A,„ I Vnfr) 



2 8 n(r 

1 /■ ,9 /■ ,9 ,n(r)n(r' 




+ ^Jd"^J d'^'j^^ + / dhvUrMr) . (1) 

The first term in (|l|) is the Thomas- Fermi kinetic energy, the second term is a von Weizsacker- 
like correction to the kinetic energy,!! and the third term is the Dirac local exchange energy. 
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For simplicity, we neglect any correlation contribution. Following our earlier work,0 we 
choose a value of = 0.25 for the von Weizsacker coefficient. Our results are not strongly 
dependent on this value since the systems studied are primarily in the Thomas- Fermi regime. 
A variation of Eq. (0) with respect to the density leads to the Euler-Lagrange equation 

5E\n] 

where the Lagrange multiplier /z (chemical potential) serves to fix the total number of elec- 



trons, A^. Introducing the so-called von Weiszacker wave function, iP{y) = yn{r), into 
Eq. we obtain 

- ^VV(r) + t;efr(r)^(r) = /.^(r) , (3) 
where the effective potential is given by 




Weff(r) = 7r^2(r) - \ -tp{r) + 0(r) + Ve^^) ■ (4) 



Here, (j){r) = J dr'n{r')/\r—r'\ is the electrostatic potential arising from the electronic density 
n{r), and Vext(i') is the externally imposed potential. The latter is assumed to be periodic 
with periodicities and ay in the x and y directions, respectively, and to have inversion 
symmetry about the center of the unit cell. 

Eqs. (|^) and (^) are reminiscent of the usual Kohn-Sham equations except that only 
a single electronic orbital need be calculated. It should also be noted that the effective 
potential Ves includes a term coming from the Thomas-Fermi kinetic energy which does 
not appear in the usual Kohn-Sham potential. The required solution to Eq. (|^) is the 
ground state von Weiszacker wave function, ipQ, which is obtained by imposing the boundary 
condition 



n ■ V^o = , (5) 

at the edges of the unit cell. Here, n is a unit normal vector. To complete the specification 
of ipo we must also impose the normalization 

J dh no(r) = N , (6) 

which fixes the chemical potential fi. 

The self-consistent solution of Eqs. (|^) and can be obtained by direct iteration,!! but 
care must be taken to avoid numerical instabilities associated with charge fiuctuations. Here 
we adopt an alternative method to avoid this difficulty. We use the method of imaginary 
time evolution defined by the equation 

^ = -(i/-/i)^, (7) 

where H is the TFDW Hamiltonian given by 



4 



A 



(8) 



Using a first-order approximation for tlie time derivative on tfie left hand side of Eq. we 
have 



ipit + 6t) = ij{t) + (yu - H)tP{t)6t , 



where /i is given by 



fx 



(9) 



(10) 



Once an initial guess for ip{t = 0) has been made, Eq. (|]) provides an updated ipi^t) 
which is used in the next iteration to determine an updated H and fi, and so on until 
convergence. This scheme preserves the normalization of the wave function to 0(5t^), but 
the wave function is nonetheless normalized to unity at each iteration. A combination of 
real-space and Fourier-space representations is used to evaluate the effect of H on ip. We 
have found this method of generating ipo to be very stable against the charge fluctuation 
induced Coulomb instability. 



B. Hydrodynamic Equations 

A theory for the collective excitations of a modulated 2DEG is constructed by treating 
the electronic system as an inviscid, charged "classical" fluid. In particular, we adopt the 
hydrodynamic equationslii'ii 



'dt 



V ■ (nv) = , 



(11) 



and 



n 



(9v „ 
_ + v.Vv 



where the total force acting on a fluid element, F*°* 
force 



(r,t) = -V 



Ves{r,t) 



nF'°\ (12) 
pint _|_ jpcxt^ consists of the internal 

A^V2^(r,t)' ^^^^ 



2 i;{r,t) 

and the force due to external electromagnetic fields (c = 1) 

We wish to show here that these hydrodynamic equations can be obtained within a 
Lagrangian formulation. We suppose that the external electromagnetic field is defined in 
terms of a scalar potential c/)*'^* and a vector potential A*^^* by 
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f)\ext 

gext ^ _V0Cxt gext ^ y ^ ^ext ^ (-^5^ 

and that the velocity field v(r,t) is defined by 

Y{r,t) = Vg,{r,t) + g2{r,t)Vgs{T,t) + A'^' . (16) 

Here, the Qi^s are three independent scalar functions; gi is associated with the irrotational 
part of the fiuid fiow and g2 and g^ are introduced to represent the solenoidal part. The 
dependence on the vector potential is displayed as a separate term. 
The Lagrangian density that we propose is given by 

(v^ dgi dgA ^ . . 

/: = -(^y + ^ + ^72^j-n0 +s[n] (17) 

where e[n\ is the TFDW energy density. Application of the principle of least action 

6 [ dt [dr C = , (18) 



for each of the variables {gi ,g2,g3,n} generates four equations. It is straightforward to 
show that independent variations of the g^ yield the following equations: 

Oti 

_ + V(nv) = 0, (19) 
^ + v.Vg2 = 0, (20) 

^ + vV^?3 = 0. (21) 

Eq. (|T9|) is just the expected continuity equation. 

Now, a variation of Eq. (0) with respect to n yields 

Taking the gradient of Eq. (|22|) , along with the time differentiated form of Eq. (|16D, we 
obtain 

^ + (v.V)v + E- + vx(Vxv)-^V,3 + ^V,, + vL.-^^)=0. (23) 



Using Eq. (|T6]) to evaluate the triple cross product, we have 

V X (V X v) = (v ■ Vg3)Vg2 - (v ■ \/g2)Vg3 + v x B^^* , (24) 
and substituting this result together with Eqs. (^) and (|2l|) into Eq. (^31), we finally obtain 



6 



dv 
'dt 



V ■ Vv = -V U, 



'eflf 



2 ip 



^ _ (e""* + V X B'="*) 



(25) 



Eqs. (|T9D and (^) are equivalent to the two equations (0) and ([T2|) written down at the 
beginning of this section. These are the fundamental hydrodynamic equations of the TFDW 
theory. 

Linearizing our hydrodynamic equations about small deviations from equilibrium, viz., 



n{r,t) — >• no + dn{r,t), and retaining only first-order quantities, Eqs. ([TT|) and (0) yield 

d6n 



dt 



+ V ■ (nov) = 



and 



dt 

where the fluctuating force is given by 

A 



V X B'^ 



5F"^*(r,t) 



Sves{r,t) 



2^C 



W'5^{v,t) 



2 



6ij{r,t) 



-V/ 



(26) 



(27) 



(28) 



with 



5v,s{r, t) = 27r^/'o5^(r, t) - \ -5ij{r, t) + (50(r, t) . 



IT 



(29) 



The quantity / plays the role of an effective potential for the internal force fluctuation. The 
magnetic field B''^* will be taken to be a uniform field in the 2;-direction. 

For the moment, we shall consider the normal modes of the system which are determined 
in the absence of the driving field E'^^*. Combining Eqs. (|2^) and (|27|) leads to the following 
equation for the fluctuating part of the von Weizsacker wave function {6tp = 6n/2ilJo) 



Uoiuj^ - ^2)5^ = _1^V2(^0/) + loo {^^] (M) + tUJcWo X V/) ■ Z 



(30) 



Here, Uc = eB°^^ /m*c, is the cyclotron frequency expressed in natural units. Making explicit 
use of the periodicity of the system in both spatial directions, we can write the fluctuating 
part of the von Weiszacker wave function in the form 



5ij{r,t) 



3«(q-r-a;t) 



iGr 



(31) 



where cg is a Fourier expansion coefficient, q is a 2D Bloch wavevector restricted to the first 
Brillouin zone, and G = [211171/ 2im/ay), with m, = 0, ±1, ±2, are reciprocal lattice 
vectors. The plane wave basis functions satisfy the or t honor mality condition 



A J A 



(32) 



where A = Qij^Cty IS the area of a unit cell. Expanding the product [ipof) in a similar way, 
viz., 
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2 ipo 2 



E foe'""-' , (33) 



G 



and substituting these expansions into Eq. (|30D , we obtain the nonhnear eigenvalue equation 
u;(u;^ -u;^)cg = c^E^gg'/g' +tUcE^GG'/G' , (34) 



where 



and 



G' 



1. „,2. 1 



Bgg' = (q + G) ^GG' + -^{ves -f^)[G- G'] , (35) 

/ A,,, 



Agg' = i[q X (G - G') + (G X G')] ■ z lnV^o[G - G'] . (36) 

We adopt the notation that square braces denote a Fourier-transformed variable. 

The Fourier coefficients /g are connected to the fluctuating wavefunction expansion 
coefficients cg through Eq. ([33|), viz., 



and 



/g = E[^GG' + ^GG' + ^GG' + A»5gg']cg' 
G' 

= Y.Mgg'Cg' , (37) 

G' 



M^G' = 27rz^o'[G-G'] 
M^G' = -\/f^o[G-G'] 

M-,,^4.E^£^S2S^, (38) 



^ V'(g. + G-)2 + (g, + c^y 

are the kinetic, exchange, and Hartree matrices that arise when the tpoSves term is Fourier 
transformed. The eigenvalues u of Eq. ( |5^ give the excitation frequencies of the system, 
and the corresponding eigenvectors c can be used to determine the density fluctuation of 
the collective mode. 

Eq. (0) represents an infinite dimensional matrix problem that, for practical purposes, 
must be solved on some subset of the 2D reciprocal lattice vectors. Any truncation of the 
basis set must be checked to ensure that the results for the modes of interest do not depend on 
the number of G vectors retained. If the system is weakly modulated, the coupling between 
different reciprocal lattice vectors is likewise weak, and the matrix problem is only of a 
modest size. However, if the system is strongly modulated with a large unit cell (a > 800 
nm), the number of basis functions required to adequately describe the dynamics of the 
system can become unmanageably large. This is a possible limitation of the plane-wave 
expansion technique. 
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Since our main interest is in making contact with the experimental FIR data, we shall 
only consider the q = limit of our calculations^ In this situation, our dynamical equations 
are invariant with respect to the point group (for = ay) when Uc = 0- A nonzero mag- 
netic field lowers the symmetry of the system to C4 because of the additional term involving 
the Agg' matrix in Eq. (0). In either case, we can appeal to the inherent symmetry of the 
reciprocal space to substantially reduce the number of effective basis vectors considered in 
a calculation. The way in which this is done is discussed in Appendix ^ 



C. Power Absorption 

We now consider the response of the system to a spatially uniform radiation field incident 
normally on the sample and polarized in the x-direction: E'^^*(r,t) = |i?o(e~*'^* + e*'^*)x. 
This is the most direct way of making contact with experiments which probe the collective 
excitations of the system by means of FIR spectroscopy. 

The physically relevant quantity is the time averaged power absorption, which is given 
by§ 



1 " — ^ 72„ -ind 



-EoReJ d'rf:^{r,uj) , (39) 

where j™'^(r,u;) = —nQ{r)vx{r,w) is the induced current density and the velocity is now a 
solution of 

= + SF'°' - V X B''"* . (40) 

at 

The total fluctuating force includes the external contribution ^F*^^* = — E*^^*. In addition, 
we have introduced a phenomenological relaxation rate 7 which has the effect of giving 
the excitations a finite lifetime. This implies that any frequency factors arising from the 
time-derivative of the velocity must be replaced hy cu ^ u = cu + ij. In previous studiesli^ 
it was found that the experimental transmission data was most faithfully reproduced by a 
frequency dependent relaxation rate, 7 = 70/a;. We retain the same frequency dependence 
here. 

Noting that the current density is itself a periodic function of both spatial directions, we 
can write Eq. (|39|) as 

i^)t ^ 1 17 o„ .-ind 

A 2 

where A is the sample area and 

■indr/~i , .1 ^ I j2„ „— iG r -ind/ 



EoRef:^^[G = 0,u] , (41) 



jr[G,u;] = -y rf^re-^^-jr(r,^) (42) 
is the Fourier coefficient of the induced current. The G = component is just the average 



induced current within a unit cell. A straightforward calculation paralleling Ref. [ |23| leads 
to the following expression for this component: 
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[G = 0, ^] = E G'^MGVg' 



The net effect of including an external driving field E''^* was previously shown to convert 
the nonlinear eigenvalue problem in Eq. (0), into a set of inhomogeneous equationsil 



)cg,-ujY^ Bccfc -^^cYl ^gg'/g' = be , (44) 



where 

bG = -EoiiCoG, - LJ,Gy)MG] . (45) 



The q = solution of Eq. ( pf ) along with Eq. ( p7| ) provides a complete determination of 
jmd|-Q _ Q^|^j_ ^Y^Q substitution of Eq. (^) into Eq. (|4l|) yields the final expression for the 
FIR power absorption. 



III. QUANTUM DOTS 

In Sec. y, the formalism for studying a 2D modulated electronic system was developed 
for an arbitrary periodic confining potential fext(r). In this section, we restrict our attention 
to potentials that confine the electrons in all three spatial dimensions, thereby forming an 
array of quasi-OD quantum dots. In particular, we have in mind the experiment performed 
by Demel et a/.^ll in which an array of quantum dots was prepared from modulation-doped 
AlGaAs/GaAs heterostructures. In this experiment, resonant anticrossings were observed 
in the magnetic dispersions for arrays of dots with N = 25 and N = 210 electrons per dot. 
These anticrossings were tentatively explained in terms of a nonlocal interaction associated 
with the finite compressibility of the electron gas,i'i although subsequent work pointed to the 
importance of noncircular perturbations. Here we explore the origin of these anticrossings 
in greater detail by systematically studying confining potentials of the form 

111 

Vexti't^, G) = 2^^^ ~'~ ^^""^ -per cos (46') , (46) 

where 9 is the azimuthal angle in polar coordinates. For p = 0, we have parabolic confinement 
with uoq = \/k. Setting e = with k,p ^ introduces a quartic radial perturbation, and 
finally, if all three parameters are nonzero, the last term in ( ^61) breaks the circular symmetry 



of the potential. It is understood that the potential in Eq. ( ^61) is centered in the unit cell 
and periodically extended throughout the array. Since the dots of interest are localized 
near the center of the cell, the behaviour of the potential near the cell boundaries is of no 
consequence. 

In order to make contact with previous work on this problem, it will prove useful to 
elaborate on the various parameters appearing in Eq. (^61). To begin, we set e = 0. Following 



Ref. [|T2[) we define the anharmonicity parameter rj = pR^ /2k, where Rq is the radius of the 
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dot as determined within a classical calculation which includes only the Hartree interaction 
of the electrons. This radius is given implicitly by the equation^ 

^0 = jTrfiyT' ' ^^^^ 

where Rq is the equilibrium radius for pure parabolic confinement (hence the superscript 
'p') given byil 



Rq 



It was further shown in Ref. [ |12[ that the zero-field {B = 0) frequency of the lowest dipolar 
mode is given to an excellent approximation by 




tu'o = LJo\ 1 + -V , (49) 



where cuq is the bare harmonic frequency defined above. Thus, in the absence of a noncircular 
perturbation, a given zero-field dipole frequency is determined by a family of ujq and rj values. 
These values are adjusted to reproduce the experimentally observed frequency u'q, as well 
as other features of the observed absorption spectrum. In this way, the quantities k and 
p can be determined. We observe from our calculations that the subsequent inclusion of 
the parameter e has a negligible effect on uJq{B = 0), so that the effect of this parameter 
can be investigated independently of the other two. Typical values for r] and e used in the 
calculations lie in the range [0, 1]. 

Finally, we summarize the appropriate physical parameters for the AlGaAs/GaAs het- 
erostructures under consideration: the effective Bohr radius is 103 A and the effective 
Rydberg of energy, Ry* = e^/2eaQ, is 5.4 meV. To be as consistent as possible with the 
experimental situation,0 we focus on dot arrays with N = 210 electrons per dot and a FIR 
zero field absorption frequency of uJq = 32 cm~^. 



A. Parabolic Confinement 



As mentioned above, rj = provides a potential that defines an array of parabolically 
confined dots. The nature of the collective excitations for such a system are by now well 
established in the literature.iHii'&S In Fig. |l](a) we show the calculated ground state density 
for a square array with ax = ay = 1000 nm containing well-separated dots with = 210 
electrons per dot. The radius of the dot is about 150 nm, and for the sake of clarity, 
we show only that portion of the unit cell in which the electronic density is localized. The 
radius found here on the basis of the TFDW functional is similar to the value obtained using 
Eq. (|i8|), viz., Rq = 160 nm, with c^o = 32 cm~^ and A^ = 210. The similarity of these values 
indicates that the Hartree interaction is in fact dominating the equilibrium distribution and 
that the exchange and kinetic energies are playing a relatively minor role. The shape of 
the distribution is approximately semi-circular as found in the classical calculationll3, but 
with a smooth variation at the edge of the dot as a result of the von Weiszacker correction. 



11 



This term in the energy functional allows us to simulate the behaviour of the charge density 
expected within a fully quantum-mechanical calculation. 

It is perhaps worth emphasizing that the equilibrium calculations described above effec- 
tively correspond to isolated quantum dots. The Coulomb interactions between equilibrium 
charge densities in different cells of the array are eliminated, so that only intra-cell Coulomb 
interactions are retained in the calculation. Of course the potential within a given cell 
arises from all other charges in the system, and one should view the parabolic confining 
potential as the net effect of all these other charges. However, when calculating the magne- 
toplasmons, the inter-cell Coulomb interactions arising from the density fluctuations are not 
neglected. These multipolar interactions are required in order to simulate the response of an 
array of coupled dots. We shall examine shortly the importance of these inter-dot Coulomb 
interactions. 

In Fig. H, we show the calculated magnetic dispersion for the array of dots described 
above. To avoid cluttering the diagram, we have chosen to omit the large number of edge 
modes which start from = at -B = and lie below u = Uc- These modes make 
no observable contribution to the power absorption and are therefore of little interest in 
the present paper. The mode frequencies found here are very similar to those found in 
the classical approximationi and can be classified according to the local 5*0(2) symmetry of 
circularly symmetric confining potentials. As discussed in Appendix ^ the Lie group 5*0(2) 
has m one-dimensional irreducible representations given by x^i^P) = e*'"'^, with < < 2tt, 
and m & Z. The different branches seen crossing in Fig. ^ correspond to different irreducible 
representations of 5*0(2) (see Table |ITl|in Appendix and thus one would not expect any 
mode repulsion in the absence of some perturbation. Since the only possible perturbation in 
the present situation is the dynamic inter-cell Coulomb interaction, the results shown here 
confirm that this interaction is very weak. As we shall see, the effect of interdot interactions 
can be made apparent by reducing the interdot separation. 

Within the classical theory, a subset of the dot modes has a dispersion given by the 
expressionili'0 

a;,,± = + {u;j2y ± tuj2 . (50) 

The solid curves in Fig. |^ show some fits to Eq. (0), and demonstrate that the field depen- 
dence is well reproduced even though the TFDW values of Qj are slightly different from the 
analytic classical values. The lowest pair of curves corresponds to the two circular polariza- 
tions of the center-of-mass (CM) mode. According to the generalized Kohn theorem,!! the 
exact separability of the A^-body Hamiltonian into CM and relative coordinates for parabolic 
confinement ensures that Qi = Uq {uq = 32 cm~^ in this case), and that only these dipole 
modes couple to a uniform external electric field. This fact is illustrated in the inset to Fig. |^ 
where we show the calculated FIR power absorption for the dot array for a range of magnetic 
fields consistent with the experiment in Ref. [ The single peak at uj{B = 0) = 32 cm~^ 
reflects the fact that our model satisfies the generalized Kohn theorem. It is also clear from 
this figure that the two-peak structure appearing for B 0, corresponds to the CM-mode 
dispersion seen in Fig. ^. 

We next consider the effects of interdot Coulomb coupling on the magnetoplasmon ex- 
citations. We have already indicated that this interaction is weak for a lattice constant of 
a = 1000 nm, and we now support this assertion by observing the effect of reducing the 
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lattice constant to a = 600 nm. Fig. ^ shows a small section of the magnetic dispersion 
for this case. Although this dispersion looks similar to Fig. ^ there are some important 
differences. 

First, we notice a slight red-shift of the uJo{B = 0) frequency while the higher- lying modes 
have moved slightly in the opposite direction. The softening of the ujo{B = 0) frequency can 
be explained by invoking the interaction between the induced dipoles on each lattice site. 
For the square lattice being considered, the electric field experienced by the m-th dipole due 
to all other dipoles is given byll 



E„ 



1 



a3f- (/2 + n2)3/2 



3f 



1 



Pm+h 



(51) 



where we have assumed, consistent with Sec. 11 Q , that the polarization is in the x-direction. 
The quantities / and n identify positions on the 2D lattice, viz., r = /ax + nay, and the 
primed summation indicates that the = / = point is excluded. In the FIR regime, each 
of the dipole moments Pm has a common magnitude po? and the electric field can be written 

as 



Po 



S 



(52) 



where 



9- ^ 



31' 



- 1 



(53) 



The quantities Pm and E^ are connected through the relation pm = a{uj)Em where a{u) is 
the dipole polarizability of the dot. For parabolic confinement, the polarizability is given by 



a{uj) = Ne'^/m*{ujQ 
the CM mode: 



Using this result in Eq. (|52|) , we find the following frequency for 



1 ^2 

^0 - ^^P ' 



(54) 



where Up = AnNe'^ /em*a^ is the effective plasmon frequency. This simple result shows that 
the interdot dipole interactions decrease the absorption frequency from that of a single dot. 

Next, we note the appearance of small anticrossings in the magnetic dispersion that was 
not discernible for the well separated dots. These anticrossings are still difficult to resolve, 
so in the figure inset, we show a close-up view of the boxed region which clearly illustrates 
the mode repulsion. These anticrossings arise because the interdot coupling of the dots on 
the square lattice breaks the local circular symmetry of the confining potential. Specifically, 
the square symmetry of the lattice is lower than the local 5*0(2) symmetry of the dots. 
As shown in Appendix ^ the 1"*" and 3~ modes now transform under the same irreducible 
representation of C4 (see Table ^ in Appendix which means that the modes will repel 
(anticross) each other under a reduction of the symmetry. We stress however that these 
anticrossings are far too weak to account for the mode anticrossings observed in the Demel 
et al. experiments. 



13 



B. Radial Anharmonic Confinement 

The inclusion of a perturbative term of the form Svcxt{i^) = Eq. (^B|) preserves the 

circular symmetry of the confining potential, but allows for the mixing of different dipole 
modes with a redistribution of the FIR oscillator strength. The ground state density is 
illustrated in Fig. |l|(b) and is qualitatively similar to that found in the classical modehlli 
This density was generated with i^o = 20cm~^ and r] = 1. The circular ridge is a result 
of the perturbation. The calculated magnetic dispersion and FIR power absorption are 
shown in Fig. 

Comparing Figs. ^ and ^, we see that many of the features of the magnetic dispersion 
in the case of parabolic confinement are still present for the anharmonic potential. To 
emphasize this point, we have plotted the dispersions described by Eq. (^Dj) using VLj as an 
adjustable parameter. We see that Eq. ( pOD provides a good fit to the data, and indicates that 
the functional form of the magnetic dispersion is not sensitive to the form of the confining 
potential. In particular, a fit of this dispersion to the lowest CM modes does not in general 
provide a direct measurement of the harmonic term in the confining potential. The relative 
spacing between the CM modes and the higher modes does depend on the anharmonicity, 
but the latter are not usually observable in FIR experiments. 

The main observable difference between parabolic and radially anharmonic confinement 
is revealed by the theoretical power absorption in the inset to Fig. ^ The curves were gen- 
erated for a range of magnetic fields corresponding to the experimental situation in Ref . [ Q . 
The dominant dipolar peaks found for parabolic confinement are now accompanied by a 
weak satellite which tracks along the uj^ CM-like mode. The oscillator strength of the 
satellite is directly controlled by the parameter 77. A value of approximately 1 reproduces 
the observed oscillator strength at high fields, while the harmonic confining frequency was 
adjusted to ujq = 20 cm~^ in order to yield a. B = dipole frequency of 32 cm^^. These 
parameters imply a strong anharmonicity. The B = peak at a; = 32 cm~^ in Fig. |^ can 
therefore no longer be identified as a CM mode in the sense of parabolic confinement. The 
very appearance of the satellite is an indication that the generalized Kohn theorem is no 
longer valid. The perturbation has the effect of coupling the dipolar modes found in the 
parabolic limit Jll thereby making other dipole modes FIR active. A comparison of our cal- 
culated FIR power absorption to the experimental datai leaves little doubt that the satellite 
structure observed in the experiment has its origins in the anharmonicity of the confining 
potential. This conclusion confirms that reached by other workers using different theoretical 
approaches.§00 The power absorption thus provides a direct probe for determining the 
geometry of the confining potential. 

Finally we mention that calculations were also carried out for the reduced lattice spacing 
of 600 nm. As in the case of parabolic confinement, anticrossings in the magnetic dispersion 
arise, but the strength of these anticrossings is again too weak to account for those observed 
experimentally. We can therefore rule out interdot Coulomb interactions as a significant 
mechanism at the experimental interdot separations. 
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C. Noncircular Confinement 

We now consider an explicit symmetry breaking perturbation of the form given in Eq. ( ^Bf ) 
with 6 0. It is clear that such a perturbation is consistent with the geometry of the dots 
studied in Ref . [ |^ . The equilibrium density for the noncircular confining potential is shown 
in Fig. |l|(c). This density profile was generated with = 210, uq = 20 cm~^, r] = 1, e = 0.4, 
and ttj; = Qy = 1000 nm. The maximal radial extent of this dot is about 155 nm. It is 
clear from Fig. |l](c) that the deviation from circular symmetry is very pronounced, and one 
should not expect the same mode dispersions as obtained for arrays of circularly symmetric 
dots. 

Turning to the magnetic dispersion in Fig. ^ we do indeed find a very different in- 
field dependence of the collective modes. Specifically, we note the appearance of strong 
anticrossings between the lowest-lying modes, and progressively weaker anticrossings for 
the higher modes. These anticrossing are a consequence of the Sve^t{.i^,0) = ^per'^cos{4:9) 
perturbation in Eq. (^61), and occur when the angular momentum of the symmetric modes 
differ by an integral multiple of four. The lowest anticrossing is between a pair of m = 1 and 
m = —3 modes and is the most important since it is the structure observed in the Demel et 
al. experiment. A comparison of Fig. ^ with Fig. 2 of Ref. [ ^ shows that our results are 
in good agreement with regard to the size and location of the anticrossing. Once again, the 
solid curves in Fig. |^ are fits to Eq. (|50D and here, they serve to emphasize the difference 
between circular and noncircular confinement geometries. 

The magnitude of the gap occurring at each anticrossing is directly related to the pa- 
rameter £. The value e = 0.4 was chosen to best fit the observed transmission data@ and the 
inset to Fig. ^ shows our calculated power absorption. At low fields there is weak structure 
to the high field side of the main resonance which is probably too weak to be resolved experi- 
mentally. However, as one enters the field range of the anticrossing between 1 and 2 T, a new 
peak appears and leads to an absorption spectrum consisting of three peaks. This structure 
is most distinct at a field of 5 = 1.5 T which is in the middle of the anticrossing region. By 
B = 2.4T, the central peak has already lost most of its oscillator strength, and the power 
absorption takes on the characteristics of the power spectrum of Fig. ^, including the 
high-field satellite. This overall behaviour is entirely consistent with experiment.! 

IV. CONCLUSIONS 

In this paper, we have presented a generalization of our previous workil on TFDW hydro- 
dynamics in laterally modulated electronic systems. As a specific application, we have ex- 
amined the magnetoplasmon excitations in arrays of both circular and noncircular dots. By 
considering potentials which include both radial and noncircular anharmonic perturbations, 
we are able to fully explore the effects of geometric confinement on the magnetoplasmon 
excitations of the dots. 

In the case of circular dots, our results indicate that an radial perturbation can 
account for the satellite peak structure observed in the experiment s.Eli'lll We have also seen 
that the interdot coupling is too weak at the experimental lattice constant to give rise to 
any discernible anticrossing effects, in agreement with earlier results.EH These effects only 
become apparent when the lattice spacing is reduced to the order of the size of a dot. 
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On the other hand, the addition of a noncircular perturbation of the form r^cos(4^^) was 
shown to induce anticrossings in the magnetoplasmon dispersion that are entirely consistent 
with the experimentally observed transmission data. Specifically, the location for the onset 
of the first optically observable anticrossing, along with higher 5-field peak structure, was 
found to be in good agreement with experiment. Based on this result, we conclude that only 
an explicit noncircular symmetry in the confining potential can fully account for the mode 
anticrossings observed in the experiment s.0i 

The work presented here is applicable to a wide variety of 2D geometrical confinements. 
Our current interest lies in the application of this formalism to antidot arrays which are 
complementary structures to quantum dot arrays. Recent experimental work on these sys- 
tems has revealed a collective excitation spectrum very different from the dot arrays.Sil A 
detailed discussion of these systems will be presented elsewhere. 
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APPENDIX A: DERIVATION OF THE SYMMETRY-REDUCED DYNAMICAL 

EQUATIONS 

The square primitive cell of our problem results in a reciprocal lattice with 4-fold sym- 
metry, as illustrated in Fig. ^. The point group of the square lattice is 6*4^,, and the full 
point symmetry is €4^ ^{E,T), where {E,T) is the time-reversal symmetry group and E 
is the identity element .e2I In the context of our problem, we recall that we are primarily 
interested in the calculation of the FIR response of the system, which corresponds to the 
q = limit of the general set of inhomogeneous equations, Eq. (P^. In the absence of a 
magnetic field, these equations are invariant under ® {E,T), however the inclusion of a 
magnetic field breaks time-reversal symmetry and lowers the symmetry of the system from 
C^y C4. This is easily seen if one observes the symmetry group of the Agg' niatrix in 
Eq. (pGl). It follows that C4 is the relevant group symmetry to consider in the most general 
case (i.e., magnetoplasma excitations). 

Let us step back for a moment, and consider the full symmetry group of the square lattice, 
namely, C^y = {E, C4, C|, C|, m^, rriy, au, Cv}, where m are reflections about the x or y axes, 
and (J is a reflection about a diagonal of the square. We can construct a representation, F 
(we use this notation to remind us that we are at the F point of the first BZ), of this group 
by considering a function space spanned by the set of functions {|Gj)}. One natural choice 
would be to consider the set of functions defined by |Gj) = a~^e^^^'^ with 1 < z < 4 as shown 
in Fig. ^. If we choose as our canonical function |Gi), the other functions are obtained by 
acting with the C4 operator on the function |Gi): C4 |Gi) = |Gi+„), with IG5) = |Gi). 
Since the generators of the group C^^ are {C4, nix}, we need only obtain the representations 
for the generators of C^^ to construct the entire representation of the group. It is easy to 
show that generators of C^^ in this representation take the form 
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(Al) 



and that the characters for C^v are given by 

C41; E 2C4 



2m 2a 



1 4 



The character table for the irreducible representation of is shown in Table |. Since 
there is no four- dimensional irreducible representation of C^^, T must be reducible. It is a 
simple matter of class-wise character addition to determine that the F representation must 
be decomposed according to F = Fi © r4 © F5. 

Now, let us consider the dipole operator fj, = er. This operator generates the represen- 
tation F^ = Fi ©F5 under 6*4^,. The Fi representation is derived from the z-component of r, 
whereas the F5 representation is induced by the (x, y) components. Owing to the fact that 
we are only considering radiation polarized in the plane of the 2DEG, we can immediately 
see that the invariant subspace of F5 contains all of the dipole active modes of our problem. 
Thus, rather than concerning ourselves with solving the generalized set of equations, viz., 
Eq. (1^), we can use the symmetry of the system to solve for only those eigenvalues that 
will be of interest in a FIR mode calculation. 

The problem of projecting into the F5 representation is really a matter of block- 
diagonalizing our system of equations and picking out that block associated with F5. This 
amounts to finding a unitary transformation that will block-diagonalize each of the F matri- 
ces of the group 6*4^,. Motivated by the knowledge that under a nonzero magnetic field the 
group symmetry is C4, we consider the eigenvalues and eigenvectors of the group element 
C4 G C4. Indeed, since C4 is the generator of the point group C4, it is the only group element 
that we need to consider. A simple calculation reveals that the eigenvalues and normalized 
eigenvectors of F(C4) are 



A(3) 



A(2) 



A(4) 



—I 



1 

-1 
1 

-1 

1 

— i 
-1 
i 



(A2) 



The four normalized eigenvectors 
matrix f/, and its inverse 



U=- 



,{A) 



) are applied in this order to produce the unitary 
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For example, applying U to the group element C4 yields 



U-^CaU 



' 1 








0" 





-1 
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—i 



(A4) 



The unitary transformation matrix, U, will block-diagonalize all of the F matrices for the 
group C^y. In other words, the eigenvectors {\v)} of C4 define a symmetrized basis for F 
of C4t,. For the general case, Uc ^ 0, the symmetry is C4, and all the matrices in the F 
representation of C4 will consist of 1 x 1 blocks (i.e., fully diagonalized) . It is now clear that 
the eigenvalues A^^-* and A'-^^ are associated with the Fi and F4 irreducible representations 
of respectively. The eigenvalues A*^'^^ and A'^''^ are associated with F5; the vector space 
spanned by their eigenvectors {\v^^^), \v^^^)) is the two-dimensional irreducible representation 
of F5. 

It is of interest to know how the irreducible representations of C^v relate to those of C4. 
This can be done by appealing to the character tables for both the C^^ and C4 point groups 
(see Tables | and H). By inspection, we can see that 

Fi^Di 

F2^ D2 
F4-^ D2 

Fs^ 1^3 © ^4 . (A5) 



From Eq. ( [A5|) , we see that the two-dimensional irreducible representation F5 splits when 
the symmetry is lowered from C^y C4. Note also that and D4 are complex conjugate 
representations. The relevance of this fact is that in the absence of a magnetic field, these 
two representations are degenerate (time- reversal symmetry). The application of a magnetic 
field will lift this degeneracy, and one should expect two modes appearing from the F5 
contribution; one from each of the D's appearing in its decomposition. In the context of our 
model, the absence of a magnetic field implies that we only require one of either A^^-* or A'^^-' 
to obtain all of the FIR-active modes of the system. 

It is also of interest to determine what the relations are between the irreducible repre- 
sentations of C4 and 5*0 (2). The (compact) Lie group SO (2) is the symmetry group of the 
circular symmetric dots. If the dots are well isolated, they are unaware of the lattice, and 
the modes of the system can be classified according to their transformation under S0{2). If 
the dots are closer together (close enough to become aware of neighbouring dots), then the 
square symmetry of the lattice will break the local 5*0(2) symmetry. In Table |T|, we show 
the character table for some of the irreducible representation of 5*0(2) for the point group 
O4. What is of note is the compatibility relations between the group 50(2) and 0*4, viz., 

1-^ 
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3-^ D-i (A6) 

Notice that if we are restricting our calculations to the invariant subspace of Fs, we will 
only be projecting out the modes with odd m; the even modes are generated by projecting 
into Fi and F4. From these compatibility relations, we see that the modes 1"*" and 3~ 
will mix (anticross) under the symmetry lowering from 5*0(2) — > C4 because they both 
transform under the same irreducible representation of C4. An examination of Table [m| 
reveals that when the symmetry is strictly S0{2), these modes belong to different irreducible 
representations, and will therefore cross. The same argument can be used if there is an 
explicit symmetry breaking via the confining potential. Specifically, even if the dots are 
well separated, a confining potential that does not transform under 5*0(2) will lower the 
symmetry of the system and cause mode mixing (anticrossing) to occur. 

So far, our discussion has been restricted to an abstract square lattice. How do we apply 
these ideas to our problem? We note that our entire reciprocal lattice is made up of sets 
of reciprocal lattice vectors {IGio)}, (a = 1,...,4), which are generated by the symmetry 
operation O4 acting on |Gji)- It is therefore sufficient to only consider the reciprocal lattice 
vectors within the first quadrant of the total reciprocal space. The index i = l,...,nmax 
labels these vectors, and the set { | Gia) } will be referred to as a shell. The total number of 
shells retained in the calculation is rimax- 

The Fourier expansion of an arbitrary function /(r) can then be written as 

\f) = J2f^c.\Gic.) , (A7) 

ia 

where the amplitudes fia are the elements of a ^rimax x 1 column vector. Each set of 
amplitudes fia (a = 1, ..,4) can be expressed in terms of the eigenvectors of F(04) as 

/.. = E/fM'^ (A8) 

A 

and the orthonormality of the eigenvectors implies the inverse relation 

a 

The star indicates complex conjugation. Now, consider some operator O which is invariant 
under 0*4. The matrix representation of the equation 0\f) = \d) is given by 

E ^ic^^oPfoP = dia ■ (AlO) 



Making use of Eqs. (|A8| ) and (^^) we find 

E(^iiVf^ = rff\ (All) 

j 

where 
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4'^ = E^-^i'^' (A12) 

a 

Ol-^ = j:^i'^'0^a,,vf^ . (A13) 

a/3 

In obtaining Eq. (|A11| ) we have used the fact that the operator O is diagonal in the r(C4) 
basis. It is clear from Eq. (|A11 ) that each of the eigenvalues, A, defines an independent 
matrix problem. This is a direct consequence of the fact that we are working in the sym- 
metrized basis of the representation F of C4. Furthermore, suppose O = VQ with V and Q 
both invariant under C4. It is easy to show that 

k 

which means that for a given A, we need deal only with the shell matrices rather than the 
matrices defined over all four quadrants of the reciprocal space. 

Let us now apply what we have learned to Eq. (|44D in the q = limit. It is readily 
verified that Aq^q,/, -Bgg'; and Mgg'; are invariant under the symmetry group C4. Therefore, 
Eq. ( ^41) can be immediately cast in terms of the symmetrized functions, viz., 



- u^d^ - ^ E Bl^M^^c^ - ... E = 6f ) , (A15) 

k k 

where we recall that summations are over the shell index. In obtaining (|A15| ), we have used 
the relation 

/f) = E<M'^- (AI6) 
j 

In principle, Eq. ( |A15D must be solved for each of the four eigenvalues of C4. However, as 
we mentioned above, if we are only interested in the FIR-active modes of the system, we can 
project into the invariant subspace of F5, and only have to consider at most, A*^^^ and A*^^-*. 
This can be made more transparent as follows. The reciprocal lattice vector Gia is given by 
(see Fig. D 

Gic, = IGiil (cos(0i«),sin(0i„)) , (A17) 

where 

(a — 1)7T 

= 0.,+ ^ ^ ' , (A18) 

and 

e, = cos-^{^\ . (A19) 



The inhomogeneous vector on the right-hand-side of Eq. ( |A15D is determined from Eq. ( [A12[ ). 
Using Eq. ( |A17| ) in Eq. (|i5|), we obtain the expression hia 

bia = -Eq (iuCOsicpia) - UJcSm{(j)ia)) jCji |^o [G^] , (A20) 
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which yields 






-^\G,l\MG^l]Eo{^ + oOc) e''^ 
-i\Gii\ilJo[Ga]Eo {u - u^) e"*' 



b 



(A21) 



Therefore, as expected, we only require the X = ±i eigenvalues to fully describe the FIR 
response. In particular, the determination of the power absorption is reduced to a solution 
of Eq. ( [A15| ) for A = ±i. Using these solutions in Eq. (|43|) , we finally obtain the following 
explicit expression of the power absorption: 



A 



EoRe 



UJ + UJc 



LU 



UJ, 



J2\G,,\MGa]Y.Ml 



(\=i) (X=i) i0. 



Y.\G^,\MG^l]T.M. 



(A= 



(A22) 



To obtain the normal mode frequencies, we use the homogeneous version of Eq. ( |A15 ) 
with 7 = 0. In this case, all positive eigenvalues for X = ±i can be obtained from the full 
set of eigenvalues for X = i. The negative eigenvalues simply correspond to the positive 
eigenvalues for A = —i. This is the method used to generate the magnetic dispersion of the 
FIR-active modes shown in Figs. 2, 4 and 5. 
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FIGURES 



FIG. 1. Equilibrium density distributions for various confinement geometries: (a) parabolic 
confinement, (b) radial anharmonicity and (c) noncircular symmetry. Parameters defining the 
confinement potential are given in the text. In all cases the dots are arranged on a square lattice 
with ttx = CLy = 1000 nm, and contain = 210 electrons. 

FIG. 2. The magnetic dispersion for an array of well separated parabolically confined dots. 
The solid circles are the numerical solutions to Eq. (p^) and the solid curves are fits to Eq. (|50|). 
The figure inset shows the calculated FIR power absorption for a range of magnetic fields. The 
parameters defining the dots are the same as used to generate Fig. ||(a). 

FIG. 3. The magnetic dispersion for an array of closely spaced dots. The only difference 
between this figure and Fig. | is a reduction of the lattice periodicity: Ux = ay = 600 nm. The 
consequences of this shorter periodicity are a softening of the ojq frequency and the appearance of 
small anticrossings. Inset: A expanded view of the boxed region ~ (1.4T, 39 cm~^) clearly showing 
the mode anticrossings. 

FIG. 4. The magnetic dispersion for an array of radially anharmonic dots. The solid circles are 
the numerical solutions to Eq. (|3^ and the solid curves are fits to Eq. (|5^) . The parameters are the 
same as those used to generate Fig. |l](b). The inset shows the calculated FIR power absorption. 

FIG. 5. The magnetic dispersion for an array of noncircular dots. The solid circles are the 
numerical solutions to Eq.(p^) and the solid curves are fits to Eq. (|50D. The same parameters as in 
Fig. |l](c) are used to generate this figure. The inset shows the calculated FIR power absorption. The 
size and position of the pronounced anticrossing at w ~ 40 cm^ (|) are in quantitative agreement 
with experiment. For B > 3T, the power absorption is almost indistinguishable from Fig. ^ 

FIG. 6. A schematic representation of the truncated 2D reciprocal lattice with point group 
symmetry C^y. The arrows are an example of a shell of vectors {|Gq)} which are generated from 
Gi by the symmetry operation C4. The solid circles illustrate the set of G vectors used in the 
symmetry-reduced calculation; a square array of points is retained to facilitate the use of fast 
Fourier transforms. 
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TABLES 





E 2C4 


CI 




2m 2a 




1 1 






1 1 


r2 


1 -1 






1 -1 




1 1 






-1 -1 


r4 


1 -1 






-1 1 




2 









rp\ OT T 


Character table lor the irreducible 1 


"epresentations of the point group 


C4 


E C4 




c! 


c| 




1 1 




1 


1 


D2 


1 -1 




1 


-1 


D3 


1 z 




-1 


—i 


L»4 


1 




-1 


i 


T'A RT T? TT 


Character table for the irreducible 


representations of the point group C4. 




E C4 




c! 


Ci 


■~ 1 

0+ 


1 1 




1 


1 


0~ 


1 1 




1 


1 


1+ 


1 i 




-1 


—i 


1- 


1 —i 




-1 


i 


2+ 


1 -1 




1 


-1 


2- 


1 -1 




1 


-1 


3+ 


1 —i 




-1 


i 


3- 


1 i 




-1 


—i 



TABLE III. Compound characters for some of the SO {2) representations of the point group C4. 
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This figure "fig2.jpg" is available in "jpg" format from: 
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This figure "fig4.jpg" is available in "jpg" format from: 
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This figure "fig5.jpg" is available in "jpg" format from: 
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This figure "fig6.jpg" is available in "jpg" format from: 
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